Topological versus rheological entanglement length in primitive path analysis protocols

Primitive path analysis algorithms are now routinely employed to analyze entanglements in computer simulations of polymeric systems, but different analysis protocols result in different estimates of the entanglement length, N_e. Here we argue that standard PPA measures the rheological entanglement length, typically employed by tube models and relevant to quantitative comparisons with experiment, while codes like Z or CReTA also determine the topological entanglement length. For loosely entangled systems, a simple analogy between between phantom networks and the mesh of entangled primitive paths suggests a factor of two between the two numbers. This result is in excellent agreement with reported values for poly-ethylene, poly-butadiene and bead-spring polymer melts.Comment: 3 pages, no figure

Net energy analysis of solar and conventional domestic hot water systems in Melbourne, Australia

It is commonly assumed that solar hot water systems save energy and reduce greenhouse gas emissions. Very rarely has the life-cycle energy requirements of solar hot water systems been analysed, including their embodied energy. The extent to which solar hot water systems save energy compared to conventional systems in Melbourne, Australia, is shown through a comparative net energy analysis. The solar systems provided a net energy saving compared to the conventional systems after 0.5 to 2 years, for electricity and gas systems respectively.<br /

Static and dynamic properties of large polymer melts in equilibrium

We present a detailed study of the static and dynamic behavior of long semiflexible polymer chains in a melt. Starting from previously obtained fully equilibrated high molecular weight polymer melts [{\it Zhang et al.} ACS Macro Lett. 3, 198 (2014)] we investigate their static and dynamic scaling behavior as predicted by theory. We find that for semiflexible chains in a melt, results of the mean square internal distance, the probability distributions of the end-to-end distance, and the chain structure factor are well described by theoretical predictions for ideal chains. We examine the motion of monomers and chains by molecular dynamics simulations using the ESPResSo++ package. The scaling predictions of the mean squared displacement of inner monomers, center of mass, and relations between them based on the Rouse and the reptation theory are verified, and related characteristic relaxation times are determined. Finally we give evidence that the entanglement length $N_{e,PPA}$ as determined by a primitive path analysis (PPA) predicts a plateau modulus, $G_N^0=\frac{4}{5}(\rho k_BT/N_e)$, consistent with stresses obtained from the Green-Kubo relation. These comprehensively characterized equilibrium structures, which offer a good compromise between flexibility, small $N_e$, computational efficiency, and small deviations from ideality provide ideal starting states for future non-equilibrium studies.Comment: 13 pages, 10 figures, to be published in J. Chem. Phys. (2016

Smectic-C tilt under shear in Smectic-A elastomers

Stenull and Lubensky [Phys. Rev. E {\bf 76}, 011706 (2007)] have argued that shear strain and tilt of the director relative to the layer normal are coupled in smectic elastomers and that the imposition of one necessarily leads to the development of the other. This means, in particular, that a Smectic-A elastomer subjected to a simple shear will develop Smectic-C-like tilt of the director. Recently, Kramer and Finkelmann [arXiv:0708.2024, Phys. Rev. E {\bf 78}, 021704 (2008)] performed shear experiments on Smectic-A elastomers using two different shear geometries. One of the experiments, which implements simple shear, produces clear evidence for the development of Smectic-C-like tilt. Here, we generalize a model for smectic elastomers introduced by Adams and Warner [Phys. Rev. E {\bf 71}, 021708 (2005)] and use it to study the magnitude of Smectic-C-like tilt under shear for the two geometries investigated by Kramer and Finkelmann. Using reasonable estimates of model parameters, we estimate the tilt angle for both geometries, and we compare our estimates to the experimental results. The other shear geometry is problematic since it introduces additional in-plane compressions in a sheet-like sample, thus inducing instabilities that we discuss.Comment: 8 pages, 5 figure

Stretching Semiflexible Polymer Chains: Evidence for the Importance of Excluded Volume Effects from Monte Carlo Simulation

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice ($d=3$ dimensions) and square lattice ($d=2$ dimensions), varying chain stiffness by an energy penalty $\epsilon_b$ for chain bending. In the absence of excluded volume interactions, the persistence length $\ell_p$ of the polymers would then simply be $\ell_p=\ell_b(2d-2)^{-1}q_b^{-1}$ with $q_b= \exp(-\epsilon_b/k_BT)$, the bond length $\ell_b$ being the lattice spacing, and $k_BT$ is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both $q_b$ and the chain length $N$ are varied over a wide range $(0.005 \leq q_b \leq 1, \; N \leq 50000$), and also a stretching force $f$ is applied to one chain end (fixing the other end at the origin). In the absence of this force, in $d=2$ a single crossover from rod-like behavior (for contour lengths less than $\ell_p$) to swollen coils occurs, invalidating the Kratky-Porod model, while in $d=3$ a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in $d=2$, while theories based on the Kratky-Porod model are found to work in $d=3$ for stiff chains in an intermediate regime of chain extensions. While for $q_b \ll 1$ in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g. bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.Comment: 23 pages, 17 figures, 2 tables, to be published in J. Chem. Phys. (2011

Optimized cross-slot flow geometry for microfluidic extension rheometry

A precision-machined cross-slot flow geometry with a shape that has been optimized by numerical simulation of the fluid kinematics is fabricated and used to measure the extensional viscosity of a dilute polymer solution. Full-field birefringence microscopy is used to monitor the evolution and growth of macromolecular anisotropy along the stagnation point streamline, and we observe the formation of a strong and uniform birefringent strand when the dimensionless flow strength exceeds a critical Weissenberg number Wicrit 0:5. Birefringence and bulk pressure drop measurements provide self consistent estimates of the planar extensional viscosity of the fluid over a wide range of deformation rates (26 s1 "_ 435 s1) and are also in close agreement with numerical simulations performed by using a finitely extensible nonlinear elastic dumbbell model

Oscillating Fracture in Rubber

We have found an oscillating instability of fast-running cracks in thin rubber sheets. A well-defined transition from straight to oscillating cracks occurs as the amount of biaxial strain increases. Measurements of the amplitude and wavelength of the oscillation near the onset of this instability indicate that the instability is a Hopf bifurcation

Deforming glassy polystyrene: Influence of pressure, thermal history, and deformation mode on yielding and hardening

The toughness of a polymer glass is determined by the interplay of yielding, strain softening, and strain hardening. Molecular-dynamics simulations of a typical polymer glass, atactic polystyrene, under the influence of active deformation have been carried out to enlighten these processes. It is observed that the dominant interaction for the yield peak is of interchain nature and for the strain hardening of intrachain nature. A connection is made with the microscopic cage-to-cage motion. It is found that the deformation does not lead to complete erasure of the thermal history but that differences persist at large length scales. Also we find that the strain-hardening modulus increases with increasing external pressure. This new observation cannot be explained by current theories such as the one based on the entanglement picture and the inclusion of this effect will lead to an improvement in constitutive modeling

A Model for the Elasticity of Compressed Emulsions

We present a new model to describe the unusual elastic properties of compressed emulsions. The response of a single droplet under compression is investigated numerically for different Wigner-Seitz cells. The response is softer than harmonic, and depends on the coordination number of the droplet. Using these results, we propose a new effective inter-droplet potential which is used to determine the elastic response of a monodisperse collection of disordered droplets as a function of volume fraction. Our results are in excellent agreement with recent experiments. This suggests that anharmonicity, together with disorder, are responsible for the quasi-linear increase of $G$ and $\Pi$ observed at $\varphi_c$.Comment: RevTeX with psfig-included figures and a galley macr